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Our numerical calculations on small carbon clusters point to the existence of a metastable three- 
dimensional eight-atom cluster Cg which has a shape of a six-atom triangular prism with two excess 
atoms above and below its bases. We gave this cluster the name "prismane". The binding energy of 
the prismane equals to 5.1 eV/atom, i.e., is 0.45 eV/atom lower than the binding energy of the stable 
one-dimensional eight-atom cluster and 2.3 eV/atom lower than the binding energy of the bulk graphite 
or diamond. Molecular dynamics simulations give evidence for a rather high stability of the prismane, the 
activation energy for a prismane decay being about 0.8 eV. The prismane lifetime increases rapidly as the 
temperature decreases indicating a possibility of experimental observation of this cluster. 

Carbon is known to form a rich variety of crystal structures due to its ability to exist in different valence states. 
As a result, coordination numbers in carbon compounds range from two (e.g. carbyne), through three (e.g. graphite) 
to four (e.g. diamond), while typical ("favorable") values of angles between covalent bonds are 180° in carbyne, 120° 
within graphite layers, and 109°28' in diamond. Recent discovery of the Cgo molecule [Q and synthesis of Ceo clusters 
in macroscopic quantities Q have stimulated a renewed interest in carbon nanostructures ||] . A carbon nanostructure 
may be viewed as a graphite layer transformed into a tube (nanotube), a ball (fulleren) etc, with an exception that 
such nanostructures usually comprise pentagons along with hexagons. 

A discovery of the cubane CgHg Q appeared to be very important from a fundamental viewpoint since the carbon 
atoms in the molecule CgHg are located in the apexes of a cube, so that the angles between C-C bonds equal to 90°, 
in contrast to the majority of carbon compounds. As a consequence of "energetically unfavorable" bond angles, the 
molecule CsHg is metastable and accumulates a considerable amount of energy [Q. Its cubic structure is stabilized by 
the corner hydrogen atoms. This rises the question as to whether there exist stable or metastable three-dimensional 
clusters that are composed of carbon atoms only and have 90° bond angles. 

Another interesting question discussed in the literature J5|-|| is "What is the minimum number of atoms in stable 
and metastable three-dimensional carbon clusters?" . It must be emphasized that while there is a substantial progress 
in studies of relatively large carbon nanoparticles Cat composed of N ~ 100 carbon atoms [0, the controversy still 
remains concerning the structure and energetics of small carbon clusters with N ~ 10, e.g., whether the stable C4, C6, 
and Cg clusters are cyclic or linear || . Three-dimensional carbon clusters (cages) are believed to be stable for N > 20 
only while for clusters with TV < 20 the stable structures are one-dimensional linear chains or monocyclic rings 

A characteristic feature of small carbon clusters is the existence of metastable states whose binding energies 
are lower than the binding energy of the stable cluster with the same N. For N < 14, the metastable structures are 
either one-dimensional (chains and rings) or two-dimensional (graphite flakes), see, e.g., Ref. ||. To our knowledge, 
up to now there was no experimental and theoretical evidence for neither stable nor metastable three-dimensional 
carbon clusters with N < 14 (the occurrence of metastable cage C14 has been reported in Ref. [|| based on the results 
of density functional calculations). 

The purpose of this work was to search for metastable three-dimensional clusters Cn with N < 14 by means of 
numerical simulation. One of such clusters, the cluster Cg, was found. Here we report the numerically calculated 
' structural and energetical characteristics of this cluster whose lifetime appears to be surprisingly long due to rather 
high value of its activation energy, about 0.8 eV. 

We have carried out TBMD (tight-binding molecular dynamics) simulations of small carbon clusters making use 
of a transferable tight-binding potential recently developed for carbon by Xu et al. pH. This numerical technique 
allows one to calculate the total energy of the cluster having an arbitrary atomic configuration. It had been proven to 
reproduce accurately the energy-versus-volume diagram of carbon polytypes and to give a good description of both 
small clusters and bulk structures of carbon RIl . We have checked that this technique describes the structure and 
energetics of small carbon clusters quite well, the difference in bond lengths and binding energies between our results 
and available ab initio calculations § usually did not exceed 10%. 

The binding (cohesive) energy E co h{N) of A-atom cluster Cat has been determined as (see, e.g., [||) 

E coh (N)=NE{l)-E(N), (1) 
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where E(N) is the total energy of the cluster, E(l) is the energy of an isolated carbon atom. The positive value of 
E co h points to the stability of the cluster with respect to its fragmentation into N carbon atoms. At a given N, there 
may exist several atomic configurations having E co h > 0. The cluster with the highest value of E co h is stable while 
the clusters with lower (but positive) values of E co h are metastable. A metastable cluster can transform to the stable, 
energetically favorable configuration. A characteristic time of such a process (the lifetime r) depends on the height 
of the energy barrier separating the metastable and stable configurations. 

We have thoroughly analyzed a number of atomic configurations as possible candidates to metastable three- 
dimensional carbon structures. However, all of them (with one exception that constitutes the essence of this paper, see 
below) appeared to be unstable and transformed to one-dimensional clusters or decayed into small fragments even at 
T — » 0. In particular, we have checked for a possibility of existence of a three-dimensional cluster Cs having the cubic 
structure (by analogy with the cubane CsH 8 Q), but we have found that such a cluster is unstable. Nevertheless, we 
have discovered that the eight-atom cluster Cs can exist as a three-dimensional cluster having a shape of a six-atom 
triangular prism with two excess top atoms above and below its bases. We gave this cluster the name "prismane". It 
is shown in Fig. 1. 

The binding energy of the prismane equals to E co h{8)/8 — 5.1 eV/atom. This is 0.45 eV/atom lower than the 
binding energy of the stable one-dimensional eight-atom cluster and 2.3 eV/atom lower than the binding energy of 
the bulk graphite or diamond. Hence, the prismane Cs is metastable. The lengths of C-C bonds equal to 2.31 A 
within each base, 1.28 A between the two bases, and 1.47 A between each base and the nearest top atom, see Fig. 1. 
We note that the angles between C-C bonds within the bases and C-C bonds connecting the bases are equal to 90°. 
This is one reason for a lower value of the binding energy of the prismane relative to the binding energy of the stable 
one-dimensional eight-atom cluster. Meanwhile, the angles between C-C bonds connecting the top atoms with bases 
of the prism equal to 104°, being close to the value of bond angle in bulk diamond (109°28'). 

In order to determine the energy barrier separating the metastable prismane structure and the stable chain structure, 
we have carried out the molecular dynamics simulations of the prismane decay for different values of the initial 
temperature Tj n j. The time of one molecular dynamics step constituted to — 2.72 • 10~ 16 s, about one percent of 
oscillation period of the dimer C2 . The temperature T of the cluster was calculated after each 500 molecular dynamics 
steps according to the formula 

^k B T = (E km ), (2) 

where k B is the Boltzmann constant, (Eki n ) is the kinetic energy per atom averaged over a period of time At = 500io- 
Such a time averaging enables one to avoid strong fluctuations of T stemming from the small number of atoms in the 
cluster under investigation. The initial temperature Tj n j was determined from Eq. (||) at < t < 500io- 

On general grounds, one would expect the probability of the cluster decay in a unit of time, W, to be given by a 
statistical formula 



W = W eM-E a /k B T ini ), (3) 

where the factor Wq has dimensionality s _1 , and E a is the activation energy, i.e., the height of the energy barrier 
separating the given metastable state from the stable or some other metastable state of the cluster. The cluster 
lifetime r may be defined as 

t = 1/W = t cxp(E a /k B T lm ), (4) 

where To = 1/Wq. Molecular dynamics simulations allow one to determine the lifetime t directly at any particular 
value of Ti n i as the time of cluster decay into another atomic configuration. It is convenient to go from the cluster 
lifetime r to a critical number of molecular dynamics steps N c = r/to which it takes for the cluster to decay: 

N c = N exp(E a /k B T ini ), (5) 

where N = T /t Q . 

Figure 2 shows a typical example of the binding energy per atom, E co h{8)/8, and the temperature T of the prismane 
as functions of the number of molecular dynamics steps N for the case T ini — 950 K. One can see that E co h(8)/8 
and T fluctuate near the values 5.0 eV/atom and 950 K respectively over N m 23500 molecular dynamics steps. The 
prismane preserves its form during this period of time. At N c = 23600 ± 100 the prismane decays into the stable 
eight-atom chain. The temperature rises up to 2600 K. Strong oscillation of the atoms within the chain with respect to 
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one another brings about the variation of E co h(8)/8 in time from 4.9 eV/atom up to the binding energy of equilibrium 
chain, 5.55 eV/atom. 

Note however that the lifetime r = N c to = 6.4- 10 -12 s is rather short on a macroscopic scale, while we are interested 
in phenomena which happen on time scales of seconds or even months and years. According to Eq. (|4|), such long 
lifetimes are expected for temperatures much lower than the value = 950 K used in simulations shown in Fig. 2. 
But those lifetimes are impossible to achieve by means of direct computer simulations. 

To overcome the problem, we have calculated numerically the critical values of N c for several values of Tj„j and 
fitted the numerical data by Eq. (||). The results are shown in Fig. 3. We stress that the process of cluster decay is 
probabilistic in nature. Hence, the lifetime r is not uniquely determined at a given value of Ti„i. As a result, the data 
for ln(iV c ) as a function of 1/Ti n i are somewhat scattered, see Fig. 3. Nevertheless one can see that the numerical 
data can be well fitted by a straight line over rather wide ranges of N c = 300 -j- 206000 and Ti„i = 890 -r- 1700 K values, 
in accordance with Eq. (||). The coefficients of this fit give the values of iVo = 2.1 and E a /ks = 9500 K in Eq. (Q). 

We note that the activation energy E a = 0.82 eV is very large, thus resulting in relatively high stability of the 
prismane. Indeed, since the dependence of r on Ti„i is exponentially strong, a decrease in Tj n j leads to a rapid increase 
in r, so that r ~ 10/zs at Tj„j = 400 K. An important remark is in order here. Since the total energy of the cluster is 
conserved in our simulations, the cluster is unable to decay if the maximum attainable (at a given value of Ti„i) kinetic 
energy E^ x of the cluster is less than E a . In its turn, E^ x is twice the time- averaged kinetic energy at the initial 
stage of cluster evolution, i.e., E™°^ = 2 • 8 • \k B T im = 24k B T ini , see Eq. (|). Hence, at T lnl < E a /24k B « 400K the 
lifetime of the cluster equals to infinity. Thus, in general, an extrapolation of r(Ti n i) curve to the range Ti„i < 400K 
is incorrect. 

The lifetime r is extremely sensitive to the value of the activation energy E a since the latter appears in the exponent, 
see Eq. (||). Hence, even a minor change in E a will result in a substantial variation of r at a given Tj„j. In order to 
refine the value of E a it is necessary to calculate the critical values of N c for a greater number of initial temperatures 
Ti n i. However, we believe it is unlikely that such a refinement will cause the value of E a to change significantly. 

Finally, it should be stressed that we have confirmed the existence of the three-dimensional metastable prismane 
structure by making use of other computer codes, e.g., MOPAC and that based on the empirical interatomic potential 
proposed for carbon systems by Tersoff JL0| . The overall shape of the cluster (Fig. 1) appeared to be the same, 
while the bond lengths and the binding energy were to some extent different. So, the prismane Cs certainly is not 
just an artifact of a specific simulation technique used in this study, while the values of the prismane binding energy, 
activation energy, and lifetimes at different temperatures may be refined by means of more sophisticated calculations. 

In conclusion, we have predicted the existence of a three-dimensional cluster Cs (prismane) which is the smallest 
three-dimensional carbon cluster found so far experimentally or theoretically. This cluster appears to be metastable 
and has the binding energy 0.45 eV/atom below the binding energy of the stable one-dimensional eight-atom cluster. 
However, molecular dynamics simulations point to a rather high stability of the prismane suggesting that this cluster 
may be observed experimentally. 

We are grateful to A. V. Krasheninnikov and N. E. L'vov for the help in preparation of the manuscript. The work 
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Figure captions 

Fig. 1. Prismanc C 8 . Bond lengths: d AB = 2.31 A, d AC = 1.28 A, d AD = 1.47 A. Bond angles: LBAC = 90°, 
LADB = 104°, LABE = 60° 

Fig. 2. Cohesive energy per atom E co h(8)/8 (closed circles) and temperature T (open circles) of the prismane Cs 
versus the number of molecular dynamics steps N. The time of one step to = 2.72 • 10~ 16 s. Initial temperature 
T ini = 950 K. The prismane decays at N c = 23500 23700. The lifetime r = N c t a = 6.4 • 10~ 12 s. 

Fig. 3. Plot of the logarithm of the critical number of molecular dynamics steps N c (corresponding to the prismane 
decay) versus the inverse initial temperature 1/Tj„j. Circles are the results of numerical calculations. The solid line 
is the least-square fit \n(N c ) = 0.7263 + 9496/T ini , where T ini is measured in K. 
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